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We develop a consistent method for calculating non-equilibrium Green's functions for a nano- 
sized dot coupled to electron reservoirs by tunneling. The leads are generally at different chemical 
potentials (non-equilibrium), and the dot may be a single molecule characterized by its own vibra- 
tional frequency u v h- We carry out a Keldysh functional analysis of the effective Green's functions 
of the local electron in the dot, and also of the vibration mode, in the Kondo regime. Finally, the 
tunneling current is calculated as a function of bias using non-equilibrium Green's function. We find 
stepwise increase of the current when bias exceeds even multiples of the phonon frequency. 



I. INTRODUCTION 

Observing electron transport through nano-scale struc- 
tures, often consisting of a single molecule, is an active 
scientific endeavour nowadays [1—5]. With the ability to 
fabricate, and to control the conditions of, such small- 
scale objects, we have the unique opportunity to explore 
the subtle many-body phenomena under non-equilibrium 
conditions. A good example of this will be the Kondo 
phenomenon occurring in a quantum dot flanked by a 
current source and drain [2-5]. Fixing the leads' chemi- 
cal potentials at different values lead to "non-equilibrium 
Kondo phenomena" . 

Unlike the solid-state quantum dot that has been stud- 
ied extensively in the past few years [4,5], a molecular 
dot contains a new degree of freedom associated with a 
conformational change. Such mode can be typically de- 
scribed by an optical phonon of frequency Lu p h- As a 
consequence, the tunneling current exhibits a step-wise 
increase as the bias voltage passes through an integer 
multiple of the vibration frequency u p h [1]. The coupling 
of phonon mode to electron transport is an interesting 
feature absent in a more conventional quantum dot. 

It is conceivable that the same electron-phonon cou- 
pling takes place in the Kondo regime, i.e. when the 
electrical transport is carried out by the Abrikosov-Suhl 
resonance state formed across the dot and the leads. A 
signature of phonon-assisted Kondo peak has been re- 
ported in a recent experiment using Cqq molecule as a 
tunneling bridge [3]. On the theoretical side, Lee and 
Choi considered the modification of the Kondo tempera- 
ture scale Tk when a phonon mode is present in the dot 
and found that, due to the decreasing value of effective 
Hubbard parameter U, one should expect increased Tk 
[6]. Diagrammatic treatment of the Kondo effect in the 
presence of a.c. perturbation and under non-equilibrium 
is also present in the literature [7]. 

Despite these early achievements, more progress can 
be made to develop easy-to-use techniques for sorting 
out non-equilibrium many-body phenomena. For ex- 
ample, the essence of Kondo resonance in equilibrium 



can be obtained straightforwardly in the functional in- 
tegral formalism developed by Read and Newns (RN) 
[8] . Here, a slave-boson representation of the free energy 
of the problem is expressed as an imaginary-time path- 
integral and subsequently the stationary-phase condition 
is applied to look for the mean-field state. The mean- 
field solution precisely corresponds to the Abrikosov-Suhl 
resonance state. On the other hand, the first pioneer- 
ing study of the non-equilibrium Kondo physics in the 
Anderson model taken up by Meir, Wingreen, and Lee 
used non-equilibrium generalization of the non-crossing- 
approximation (NCA) and equation-of-motion (EOM) 
techniques [9] . Aguado and Langreth recently introduced 
the non-equilibrium generalization of the slave-boson the- 
ory of the Kondo resonance in their study of the transport 
through double quantum dot while avoiding the direct 
use of the free-energy minimization principle of RN [10]. 

In this paper, we take the first step towards a con- 
sistent generalization of the functional integral theory of 
RN to the case of non-equilibrium Kondo physics. In- 
stead of the imaginary time integral, we write down the 
real-time action over the Keldysh contour. Integrating 
out the lead fermion degrees of freedom, the effective ac- 
tion for the local fermion in the dot is derived that agree 
with the expressions obtained by earlier authors [9,10]. 
The method is then applied to study the phonon-couplcd 
Kondo resonance under non-equilibrium. We obtain the 
linear damping rate of the phonon mode due to coupling 
to the current at finite bias V. Furthermore, current vs. 
voltage formula is obtained for this case that exhibits 
the step- wise increase in the current, at even multiples of 
LUph, for strong electron-phonon coupling and for a nar- 
row width of the resonance spectrum. 

Before we go into the details of the result it is worth- 
while to note a similar effort by Komnik and Gogolin, 
who combined effective action approach with the Keldysh 
technique to work out the mean-field theory of the An- 
derson impurity model under non-equilibrium condition 
[11]. 

The following sections are organized as follows. Sec- 
tion II introduces the main formalism through discus- 
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sion of the non-equilibrium Kondo resonance. Section 
III is the generalization to phonon-coupled Kondo prob- 
lem, in non-equilibrium. Current-voltage expression for 
the phonon-coupled Kondo model is derived in section 
IV. We close with a summary and outlook in section V. 

II. NON-EQUILIBRIUM KONDO RESONANCE 

We consider the simplest non-equilibrium case of two 
leads with chemical potentials pl and pa, connected 
by tunneling to a quantum dot. We employ the model 
Hamiltonian H = Hl + Hr + Hd + Ht where 

H MR) = H( £fc ~ ^HK)) C L(R)ka C UR)ka 

H d = -edn d + Un d (n d - 1) 



H T = V L Y J (4 ka d, + d+ a c Lka ) +(£,—> R), (2.1) 

ka 

to study such system. We denote the local fermion oc- 
cupation by rid = Y]„ d~^d a . There is no notion of well- 
defined set of eigenstates for a non-equilibrium problem 
such as this one. Consequently, one is unable to write 
down the free energy for the problem according to the 
textbook recipe. This is why, we believe, naively writing 
down the imaginary-time action for p^ ^ p R leads to 
incorrect results. Instead, we write down the Lagrangian 
density in real-time coordinate along the Kcldysh con- 
tour [ti,tf] x [tf, U] with tf — U = T. Eventually we take 
U — ► — oo, and tf — * oo. First, the Lagrangian density 
for the Kondo problem, Eq. (2.1), in the slave-boson 
representation (U — ► oo) is written down: 

R) 

(2.2) 



^) = E C+ LkMt - <* + P L )c L ka - VL Yyik^fa + Hhc Lka \ + (L 

k k 
+ f+(id t + e d )L + b+id t b - X(f+f a + b+b - 1). 



(Through the rest of the paper we assume the spin sum- 
mation and drop X^-) This Lagrangian density is inte- 
grated along the Keldysh contour C to give the action 
S = j c dtC(t). From this we isolate the dynamics of the 
dot by integrating out the lead electrons c^ka, CRka- 

£e ff = -f+b[vl J2 Ghk + yg£ G Rk \b+u 

k k 

+ f+(id t + € d )U + b+id t b-\(f+f a +b+b-l). (2.3) 

The self-energy part is defined over the double Keldysh 
contour J c J c dtdt' while the remaining terms are over 
a single contour, J c dt. Green's functions for the lead 
electrons are given by G L ( R)k = (id t -e k + Ml(h)) -1 - At 
this stage, we take the mean-field ansatz, b(t) = b and 
X(t) — A, over the entire contour C. When this is done, 
in fact the X(b + b— 1) part drops out from the Lagrangian 
because the forward and backward integrals cancel out. 
We thus have the effective action (Gl — J2 k Glh-iGr = 
J2k GRk, A' = A - e d , V[ 



L(R) 



bV, 



>eff 



[ dtf+(t)(id t 
JC 



- f dtdt'f+(t)[V L 2 G L (t,t') 
Jc 



L(R) ) 

A')A(t) 

HL^R)]h(t')- (2.4) 



The Green's function Gum(t,t') may be labeled 

(++), (H — ), (— +)j ( ) depending on whether the time 

argument t, t' are on the upper (+) or lower (-) branch 
of the contour. In turn, we can associate each Green's 
function with G T , G < , G > and G T that are commonly 
used in the non-equilibrium theory as 



G+-(t-t') 
G~ + (t-t') 
G—(t-t') 



G<(t-t') 
G>(t-t') 

: G T (t - t') 



-%m)W)) 



In addition retarded and advanced Green's function may 
be obtained readily by G R = G T - G< , G A = G T - G> . 

We must also distinguish the fermion operators de- 
fined in the upper and lower branches separately, as 

fa+ 
fa- 

of Eq. (2.4) becomes 



F, 



In Fourier space, the effective action 



co-A'-S++ 



-cj+A'-S" 



Fa(t0), 



where we introduced the simplifying notation S ah = 
{V'*Gf + V R 2 G R b ), a, b = +, - and £ w = / ^. 

The inverse of the 2x2 matrix in the above effective 
action is the Green's function for the local fermion. Lead 
Green's functions are readily obtained 



(Gl 


Gl\ 




\G> 


Gl) 


1 = -ivp ( 



sgn(a)-,u L ) -20{hl-u) 
20(uj - p L ) sgn(w - p L ) 



, (2.5) 



assuming a constant density of states p. Similar defini- 
tions apply for the electrons in the right lead. We assume 
that the left/right leads have the same density of states. 
On inverting the matrix, we get for the local fermion 
Green's function 



G ++ (t - t') = G T (t - t 1 ) = -i(Ti>(t)W)) 



2 



G T f G< 



G> G]- ) (co-A') 2 +r 2 

w-A'-j(r L sgn(w-/i L ) + T R sga((j-fj, R )) 2i{T L 6{p, L -u) + T R 6{p, R -uj)) 

-2i(T L 9(uj-fi L ) + r R 9(uj-fi R )) i(T L sgn(uj-fi L )+T R sgn(uj-iJ, R )) 



(2.6) 



where T L = tt/^ 2 , T R = npV R and r = T L + T R . 
The non-equilibrium Green's functions obtained in this 
way are identical to those obtained by earlier authors 
[9,10]. However, we argue that the present method is 
much simpler in the sense that only the Gaussian in- 
tegration is involved in evaluating the effective Green's 
function. One should introduce two species of fermions, 
for the upper/lower branches of the contour, but can oth- 
erwise proceed with the usual theory of integration over 
quantum fields. The retarded Green's function is given 
by Gf = G T f - Gf = l/(w - A' + iT). The spectral 
function, given as the imaginary part of G*f , indicates 
the formation of a resonance state of width T, centered 
at the energy A'. 

Determination of the precise value of A' is hampered 
by the lack of precise definition of the free energy at 
hand, as discussed at the beginning of this section [12]. 
How to obtain a suitable generalization of the free energy 
minimization principle for non-equilibrium is outside the 
scope of the present paper. Instead, we resort to ref. 
[10] to conclude A' « (fiL + Hr)/2- Therefore, a slave- 
boson theory in non-equilibrium predicts a single Kondo 
resonance peak centered half-way between the chemical 
potentials of the leads as noted earlier [10]. 



III. PHONON-COUPLED KONDO RESONANCE 
AT NON-EQUILIBRIUM 

A molecular structure placed between the metallic 
leads, as well as serving as a bridge for electrons to hop 
over, may vibrate thermally, or quantum-mechanically in 
response to interaction with the passing electrons. Such 
phonon signature in the I — V characteristics has been 
first detected in the tunneling experiment with the C 60 



molecule. In this experiment, the observed phonon mode 
is that of the center-of-mass displacement of the molecule 
as a bridge [1] . A schematic drawing of the experimental 
setup is shown in Fig. 1. 

The geometrical distance between the molecule and the 
lead varies as the center-of-mass position, denoted by x, 
swings back and forth. The tunneling amplitude, being 
proportional to the overlap of the wavefunctions in the 
lead and the molecule, becomes dependent on x. The 
tunneling Hamiltonian will be accordingly modified: 



H T = V L (x) ]T c+ k J a + V R {x) V c+ 



kcr 



ka 



Rka 



d a + h.c. (3.1) 



The motion of the molecule itself is modeled with a har- 
monic oscillator of frequency co p h- In this section, we 
discuss how the phonon mode is damped due to inter- 
action with the electrons that form the resonance. The 
influence of phonons on the electronic sector has been 
considered in detail in Ref. [6]. Hereafter we adopt the 
term "phonon-Kondo model" to describe the situation 
we are interested in. 



V L (x) 




FIG. 1. Schematic setup of the phonon-coupled Kondo 
model. Left/right leads of chemical potentials fiL and hr 
are connected to a quantum dot (QD) by tunneling. The 
tunneling amplitudes Vl(r)(x) depend on the center-of-mass 
coordinate x of the QD. 

The phonon-Kondo model may be written in real-time, 



m = J2 C+ Lka(^t -e k + ^L)c Lka - bV L ( X ) YHkJo + .CcLka] + (L —> R) 



+ f+{id t -\')U + 



Wtx) 



ka 



1 2 2l 



(3.2) 



We anticipate the mean-field condition b(t) = b, X(t) = A in writing down the above Lagrangian. Getting rid of the 
lead fermions as before, the effective action becomes 



Seff = - f dtidt 2 /+(t 1 )[F/(x(t 1 ))G L (t 1 ,t 2 )F/(x(t 2 ))]/ ff (t 2 ) + (L -> R) 
Jc 

+ dtf+md t - X')U(t) + ^~J c dt[{d t xf - co 2 ph x 2 }. 



(3.3) 
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We have once again introduced the simplifying notation 
VL, R) (x) = bV L{R) (x). 

The two remaining dynamical modes in the above ef- 
fective action are the local fermions and the phonon. We 
now integrate out the local fermions to derive the effec- 
tive dynamics of the phonons, in part to demonstrate the 
utility of the Keldysh technique in describing the boson 
dynamics in non-equilibrium. We treat V' L (x) w Vq—V{x, 
V R (x) w Vq + V{x, and expand the effective action to 
quadratic order in x. In particular we have the bosonic 
part of the action 




[{d t xf - u 2 ph x 2 ] + 2iV[ 2 Tv[G f xG e x\. (3.4) 



Here G/ refers to the inverse of [id t — A' — V 2 G e ], derived 
in Eq. (2.6) of the previous section, and G e = Gl + Gr 
(See Eq. (2.5)). The self-energy portion of Eq. (3.4) can 
be written out also in the form 

2iV 2 J dtdt' 'X 1 \t) ^.--^j{t-t')X{t'), (3.5) 

where we introduced X = ( ^ + ^ , and 
^, , / D T D<\ 



K++(t) = [G T f {t)G T e (-t) + G T f {-t)G T e {t)\/2 
K+-(t) = [Gf (t)G> (-t) + G> (-t)G< (t)}/2 
K-+(t) = [G> (t)G< (-t) + G< (-t)G> (t)}/2 
K—(t) = [G T (t)G T (-t) + G T (-t)G T (t)}/2. (3.6) 

A product of Green's functions has the following iden- 
tities 

Gf (t)Gf(-t) = e(t)G>(t)G<(-t) + 6(-t)G<(t)G>(-t) 
G T 1 {t)Gl(-t) = e(t)G<(t)G>{-t) + 6{-t)G>(t)G<(-t) 

that one can use to show, in Fourier components, K ++ = 
K— = \{A + B + C + D) , K+- = \{A + B), 
K~+ = \{C + D) where A = Gf {u)G> (u> -v),B = 
£ w G>(u;)G<(u; + v),C = E u G>(o;)G<(w -v),D = 
E w G<(w)G>(w + z/). By using Eqs. (2.5) and (2.6) one 
can readily calculate A, B, G, and D and all the kernels 
in Eq. (3.6). 

On inserting the self-energy expression back into Eq. 
(3.4) and inverting the matrix, a remarkably simple ex- 
pression of the phonon Green's function emerges: 

-2iVi 2 K+- \ 
2 2 , . (3.7) 



(^) +v;\K+--K-+r 



- ^ 2 K 
-2iV[ 2 K~ + 



Retarded Green's function for non-equilibrium is the dif- 
ference 



D R = D T - D< = 



2u! p h 



9 2 



ph ^2zLo ph Vi 2 (K-+ -K+~y 

(3.8) 

using the relation 2(K+~ - K — ) = K+~ - K~+ . 
Straightforward calculation yields K h — K H = 



p I tan" 1 ^-i{l*-"-* n 



—tan 



+ P (tan- 



-tan 



-1 / fiR-v-X 



2pT 2pT 

^ L -X') 2 + r 2+ ( PR -X') 2 + r 2 , 



(3.9) 



for a pair of chemical potentials hl and p R of the 
leads and for low frequency v. Mean-field condition sets 
\fiL — A'| ~ \fJ:R — A' | ~ V/2 without the phonon coupling. 
We do not expect much modification of this conclusion 
with the phonon coupling turned on, since the strength 
of the electron-phonon coupling is significantly reduced 
to V{ = bVi <C V\ after the Kondo resonance is estab- 
lished. We then get the linear damping coefficient for the 
phonon 



1{V) 



'2 



(v/2) 2 + r 2 ' 



(3.10) 



Interestingly, we find the damping rate that is inde- 
pendent of the rcnormalization factor b 2 at zero-bias, 
j(V = 0) = (4/ir)V?/(V£ + V|). Finite V surprisingly 
reduces the damping of the phonons. 



IV. NON-EQUILIBRIUM TUNNELING 
CURRENT 



A step-wise jump in the tunneling current when mea- 
sured as a function of the bias voltage was observed in the 
experiment using the Ceo molecule as a bridge connecting 
the leads [1]. The cause of the jump is the opening up of 
an additional conducting channel as the bias voltage ex- 
ceeds integer multiples of the vibration frequency ui p h of 
the intermediate molecule. While theory of the conduc- 
tance through a vibrating molecule has been worked out 
earlier by several groups [13-15], none addressed the is- 
sue using the non-equilibrium Green's function technique 
advocated in this paper, nor in the Kondo regime as we 
do now. We show here that the pronounced features of 



4 



the experiment can indeed be shown to follow from our 
formalism. 

To this end, we consider the non-equilibrium 
Anderson-Holstein model defined by the Hamiltonian 
consisting of the left /right leads, at different chemical 
potentials pl and pr, the hybridization between the dot 
and the lead with the tunneling element Vl and Vr, now 
taken as constants, and the dot+phonon Hamiltonian 
given by H = ~e d n d + Un d (n d - 1) + uj ph b + b + gu ph (n d - 
1 ) (b + + b) . A well-known Lang-Firsov canonical transfor- 
mation of the operators, O — > e s Oe~ s , S = gn d (b + — b), 
transforms U into U — 2g 2 u; p f l , and the tunneling Hamil- 
tonian becomes 

V L ^2(4 kr rXd, + h.c.) + (L^R) (4.1) 

k 

where X = eM b - b+ ^ = e -a 2 / 2 e s b+ e a b . At the same 
time, the local electron operator assumes a composite 
form d a — > Xd a . Original electron-phonon coupling 
term has been removed in the transformed Hamiltonian. 
We further re-write the d a operator in the slave-boson 
form b + f a assuming that the effective Hubbard parame- 
ter U — 2g 2 u) p h is large enough to prohibit double occu- 
pation of the dot. 

In the Hamiltonian after the Lang-Firsov transforma- 
tion, electron-phonon coupling takes place through the 
X-operators appearing in Eq. (4.1). When the phonon 
relaxation is fast, i.e. any change in the phonon num- 
bers caused during the electron tunneling process is re- 
laxed fast, phonon number distribution will be given 
by the equilibrium value. In this case one can replace 
X w (X) equi i, « e~ 9 I 2 . The last equality holds at zero 
temperature. Now we have achieved complete separa- 
tion of the phonon and the electron dynamics in the 
transformed Hamiltonian and the electron-phonon cou- 
pling provides only the renormalization of the tunneling 
amplitudes, Vur) — * e ~ 3 ^Vur)- This is not to say 
that phonons have no effect on the dynamics of the local 
electrons. Rather, the electron Green's function, being a 
product of the phonon and the transformed local fermion 
Green's functions, will depend on the dynamics of both 
quantities. 

We proceed to calculate the Green's functions for each 
operator. The Keldysh action one obtains from Lang- 
Firsov-transformed Hamiltonian will consist of separate 
pieces for the phonon and the fermion degrees of free- 
dom according to the simplification of the previous para- 
graph. Therefore, the local electron Green's function is 
calculated as -i(d a (t)d+) = -ib 2 (X(t)X+)(f a (t)f+) = 
b 2 {X(t)X+)G>{t), where G> refers to the Green's func- 
tion of the Lang-Firsov-transformed fcrmions. The cor- 
relation function for X is obtained straightforwardly 
from the harmonic oscillator algebra: (X(t)X + ) = 
e 9 I e After Fourier transform, a simple expression 

for G> emerges: G>(uj) = b 2 e~<? £~ ±g 2n G>(uj - 



nojph)- In a similar manner we find G^(oj) = 
b 2 e-° 2 En=0 7?.9 2n Gf(u + ruv ph ). 

In calculating the current through the dot, we invoke 
the Meir-Wingreen formula [9] that requires the knowl- 
edge of the tunneling density of states p d (ui). A short ma- 
nipulation gives p d (u) = (i/2n)[G^(u)) — G^(u>)]. Cur- 
rent at zero-temperature is proportional to the integral 
of pd(u>) in the energy window pn < oj < pl, assuming 
PL — PR = V > 0. We find, up to a proportionality 
constant, the current 

2T ) + £ g -^ e (V-nu ph ) x 

{tai -. (-)"!-. ( ^)}. m 

In obtaining this expression we used the mean-field re- 
sult, A' w (p L + Pr)/2, valid in the Kondo regime. The 
current formula depends on the combined broadening pa- 
rameter T = + Tji. 

Plot of the current vs. voltage for several values of 
broadening parameter T and electron-phonon coupling 
strength g is shown in Fig. 2. It shows that the step- wise 
behavior is most pronounced for small values of T/tu p h 
and for large g. In this regard, the Kondo resonance 
regime may be an excellent venue for the observation of 
the phonon-assisted steps because of the small T value of 
the Abrikosov-Suhl state. Interestingly, the steps occur 
at even multiples of the phonon frequency, and are ab- 
sent for odd multiples of u> p h. This is to be distinguished 
from the steps at both even/odd integer multiples of uj p h 
observed in Cqq experiment. The difference is due to the 
nature of the energy level of the Kondo resonance state. 
In a resonant-tunneling model the energy level of the dot 
is fixed regardless of the bias. In the Kondo regime, how- 
ever, the resonance is located at midway between the 
chemical potentials of the leads, hence it "floats" as the 
average (pl + Pr)/2 is varied. Effectively, the relative 
energy difference between a lead and the dot level is half 
the overall bias applied. 



i a ib i c 




FIG. 2. Tunneling current vs. bias voltage (in units of u> p h) 
for several parameters (g,T/ui p h) = (a)(0.5,0.5), (b) (1,0.1) 
and (c) (1,0.01). Steps are more pronounced for a smaller 
broadening parameter F and larger electron-phonon coupling 
g, and occur at even multiples of ui ph . 

We have demonstrated that such non-equilibrium phe- 
nomena as the phonon-assisted steps in the I-V charac- 
teristics can be consistently derived using the formalism 
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developed here. 



V. SUMMARY AND OUTLOOK 

In this paper, we considered some theoretical aspects of 
transport through a prototype single-molecule device, as 
experimentally realized recently. A formalism has been 
developed that extrapolate smoothly between equilbrium 
and non-equilibrium situations. In summary, our method 
is based on writing down the real-time action over the 
Keldysh contour and integrating out the lead electrons' 
degree of freedom. This leaves us with the effective dy- 
namics of the localized electron in the dot, or for the 
phonon-coupled case, the dynamics of the oscillatory mo- 
tion of the dot. For example, we have obtained the linear 
damping rate of the phonon mode in interaction with the 
tunneling electrons in the Kondo regime. The tunneling 
current through the phonon-coupled dot has been calcu- 
lated as well; it shows step-like increase in the current vs. 
voltage behavior. 

Although the Kondo regime has been assumed through 
the discussion of this paper, it is clear that our technique 
can be applied to any out-of-cquilibrium, resonant-level 
model. This work demonstrates the utility and simplicity 
of the Keldysh functional integral approach for dealing 
with non-equilibrium transport phenomena. 
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